Search for: All records

Abstract In the past 40 years, the global annual mean surface temperature has experienced a nonuniform warming, differing from the spatially uniform warming simulated by the forced responses of large multimodel ensembles to anthropogenic forcing. Rather, it exhibits significant asymmetry between the Arctic and Antarctic, with intermittent and spatially varying warming trends along the Northern Hemisphere (NH) midlatitudes and a slight cooling in the tropical eastern Pacific. In particular, this “wavy” pattern of temperature changes over the NH midlatitudes features strong cooling over Eurasia in boreal winter. Here, we show that these nonuniform features of surface temperature changes are likely tied together by tropical eastern Pacific sea surface temperatures (SSTs), via a global atmospheric teleconnection. Using six reanalyses, we find that this teleconnection can be consistently obtained as a leading circulation mode in the past century. This tropically driven teleconnection is associated with a Pacific SST pattern resembling the interdecadal Pacific oscillation (IPO), and hereafter referred to as the IPO-related bipolar teleconnection (IPO-BT). Further, two paleo-reanalysis reconstruction datasets show that the IPO-BT is a robust recurrent mode over the past 400 and 2000 years. The IPO-BT mode may thus serve as an important internal mode that regulates high-latitude climate variability on multidecadal time scales, favoring a warming (cooling) episode in the Arctic accompanied by cooling (warming) over Eurasia and the Southern Ocean (SO). Thus, the spatial nonuniformity of recent surface temperature trends may be partially explained by the enhanced appearance of the IPO-BT mode by a transition of the IPO toward a cooling phase in the eastern Pacific in the past decades. 

In a recent work, a sensitivity analysis methodology was described that allows for a visual display of forecast sensitivity, with respect to model parameters, across a gridded forecast field. In that approach, sensitivity was assessed with respect to model parameters that are continuous in nature. Here, the analogous methodology is developed for situations involving noncontinuous (discrete or categorical) model parameters. The method is variance based, and the variances are estimated via a random-effects model based on 2^{k− p}fractional factorial designs and Graeco-Latin square designs. The development is guided by its application to model parameters in the stochastic kinetic energy backscatter scheme (SKEBS), which control perturbations at unresolved, subgrid scales. In addition to the SKEBS parameters, the effect of daily variability and replication (both, discrete factors) are also examined. The forecasts examined are for precipitation, temperature, and wind speed. In this particular application, it is found that the model parameters have a much weaker effect on the forecasts as compared to the effect of daily variability and replication, and that sensitivities, weak or strong, often have a distinctive spatial structure that reflects underlying topography and/or weather patterns. These findings caution against fine-tuning methods that disregard 1) sources of variability other than those due to model parameters, and 2) spatial structure in the forecasts.

 

Abstract. The Last Millennium Reanalysis (LMR) utilizes an ensemble methodology to assimilate paleoclimate data for the production of annually resolved climate field reconstructions of the Common Era. Two key elements are the focus of this work: the set of assimilated proxy records and the forward models that map climate variables to proxy measurements. Results based on an updated proxy database and seasonal regression-based forward models are compared to the LMR prototype, which was based on a smaller set of proxy records and simpler proxy models formulated as univariate linear regressions against annual temperature. Validation against various instrumental-era gridded analyses shows that the new reconstructions of surface air temperature and 500 hPa geopotential height are significantly improved (from 10 % to more than 100 %), while improvements in reconstruction of the Palmer Drought Severity Index are more modest. Additional experiments designed to isolate the sources of improvement reveal the importance of the updated proxy records, including coral records for improving tropical reconstructions, and tree-ring density records for temperature reconstructions, particularly in high northern latitudes. Proxy forward models that account for seasonal responses, and dependence on both temperature and moisture for tree-ring width, also contribute to improvements in reconstructed thermodynamic and hydroclimate variables in midlatitudes. The variability of temperature at multidecadal to centennial scales is also shown to be sensitive to the set of assimilated proxies, especially to the inclusion of primarily moisture-sensitive tree-ring-width records. 

As harsh weather conditions in Antarctica make it difficult to support a dense weather observing network there, it is critical to place new weather stations in locations that are optimal for a given monitoring goal. Here we demonstrate a network design algorithm that uses ensemble sensitivity to identify optimal locations for new automatic weather stations in Antarctica. We define the optimal location as one that maximizes the reduction in total variance of a given spatial field. Using WRF Model forecast output from the Antarctic Mesoscale Prediction System (AMPS), we identify the best locations for observations across the continent by considering two spatial fields: (i) the daily 0000 UTC 2-m temperature analysis field and (ii) the daily 0000 UTC 2-m air temperature 24-h forecast field. We explore the impact of spatial localization on the results, finding that a covariance length scale of 3000 km is appropriate for these metrics. We find optimal locations assuming that no stations exist on the continent (blank slate) and conditional on existing stations (CD90). In the “blank slate” scenario, the Megadunes region emerges as the most important location to both monitor temperature and reduce temperature forecast errors, with the Ronne Coast and the Siple Coast following. Results for the monitoring and forecasting metrics are similar for the CD90 subset as well, indicating that additional stations could benefit multiple performance goals. Considering the CD90 subset, Wilkes Land–Adelie Coast, Ellsworth Land, and Queen Maud Land–Interior are identified as regions to consider installing new stations for optimizing network performance.

 

Reconstructing past climates remains a difficult task because pre‐instrumental observational networks are composed of geographically sparse and noisy paleoclimate proxy records that require statistical techniques to inform complete climate fields. Traditionally, instrumental or climate model statistical relationships are used to spread information from proxy measurements to other locations and to other climate variables. Here ensembles drawn from single climate models and from combinations of multiple climate models are used to reconstruct temperature variability over the last millennium in idealized experiments. We find that reconstructions derived from multi‐model ensembles produce lower error than reconstructions from single‐model ensembles when reconstructing independent model and instrumental data. Specifically, we find the largest decreases in error over regions far from proxy locations that are often associated with large uncertainties in model physics, such as mid‐ and high‐latitude ocean and sea‐ice regions. Furthermore, we find that multi‐model ensemble reconstructions outperform single‐model reconstructions that use covariance localization. We propose that multi‐model ensembles could be used to improve paleoclimate reconstructions in time periods beyond the last millennium and for climate variables other than air temperature, such as drought metrics or sea ice variables.

 

Stochastic kinetic energy backscatter schemes (SKEBSs) are introduced in numerical weather forecast models to represent uncertainties related to unresolved subgrid‐scale processes. These schemes are formulated using a set of parameters that must be determined using physical knowledge and/or to obtain a desired outcome. Here, a methodology is developed for assessing the effect of four factors on spatial features of forecasts simulated by the SKEBS‐enabled Weather Research and Forecasting model. The four factors include two physically motivated SKEBS parameters (the determining amplitude of perturbations applied to stream function and potential temperature tendencies), a purely stochastic element (a seed used in generating random perturbations) and a factor reflecting daily variability. A simple threshold‐based approach for identifying coherent objects within forecast fields is employed, and the effect of the four factors on object features (e.g. number, size and intensity) is assessed. Four object types are examined: upper‐air jet streaks, low‐level jets, precipitation areas and frontal boundaries. The proposed method consists of a set of standard techniques in experimental design, based on the analysis of variance, tailored to sensitivity analysis. More specifically, a Latin square design is employed to reduce the number of model simulations necessary for performing the sensitivity analysis. Fixed effects and random effects models are employed to assess the main effects and the percentage of the total variability explained by the four factors. It is found that the two SKEBS parameters do not have an appreciable and/or statistically significant effect on any of the examined object features.